William I. Fine Theoretical Physics Institute 
University of Minnesota 



TPI-MINN-06/28-T 
UMN-TH-2514-06 
August 2006 



Possible new resonance at the D*D* threshold in e + e 

annihilation 

S. Dubynskiy 

School of Physics and Astronomy, University of Minnesota, 
Minneapolis, MN 55455 
and 

M.B. Voloshin 

William I. Fine Theoretical Physics Institute, University of Minnesota, 

Minneapolis, MN 55455 
and 

Institute of Theoretical and Experimental Physics, Moscow, 117218 

Abstract 

We argue that the recent CLEO-c data on e + e~ annihilation into pairs of charmed 
mesons at cm. energy around 4.0 GeV are not well described by a single resonance 
V>(4040), but can be better understood if there is an additional narrow resonance with 
mass within few MeV from the D*D* threshold. 



The strong dynamics of heavy and light quarks gives rise to rich structures near the 
thresholds of open heavy flavor states. In particular, the cross section of the e + e~ annihilation 
is well known to display an intricate behavior in the region of the thresholds for production 
of various pairs of D (D*) mesons^. The most prominent peak in the annihilation cross 
section is conventionally labeled lj as the ^(4040) resonance. It has been argued long ago|2] 
that this resonance (then called ^(4028)) might be of a 'molecular' [3. type, i.e. essentially 
made of a meson pair D*D*, on the basis of the observed great enhancement of its coupling 
to D*D* in comparison with its interaction with the other open charmed meson channels: 
DD and DD*. However it was later pointed out [4] that a similar pattern of couplings to 
meson pairs can arise for a pure charmonium state due to a kinematical cancellation in 
overlap integrals at the momentum corresponding to the meson states DD and DD*. In 
other words, the latter scheme suggests a suppression of these lighter channels rather than 
an enhancement of the coupling to D*D*. This issue has not been resolved so far due to lack 
of more detailed data on the production of the charmed mesons in e + e~ annihilation in the 
vicinity of the peak. Recently the CLEO-c experiment has produced [5J a greatly improved 
set of data from the scan of the e + e~ production cross section separately for each open 
channel with charmed meson pairs. The data display a highly nontrivial behavior in each of 
the channels in all the scanned energy range from 3970MeV to 4260 MeV, and in particular at 
energies around the ^(4040) peak. The purpose of this letter is to point out that the CLEO- 
c data in the immediate vicinity of the threshold for the D*D* meson pairs, specifically in the 
energy range from 3970 MeV to 4060 MeV, cannot be described with a reasonable statistical 
significance by a single resonance and a non-resonant background. Rather we find that the 
data suggest that in addition to a resonance with parameters close to those of ?/>(4040) there 
exists a narrow structure, with a width of not more than few MeV, at a cm. energy very 
close to 4010 -T- 4015 MeV, i.e. very close to the thresholds for the neutral vector meson pairs 
D*°D*° (4013 ± lMeV) and for the charged ones D* + D*- (4020 ± lMeV). The available 
CLEO-c data however are yet not detailed enough to fit the parameters of the resonance 
suggested by this structure. We were only able to make a very preliminary estimate of the 
e + e~ branching fraction of such resonance as B ee ~ 10~ 7 . 

If confirmed by a more detailed experimental study the suggested resonance would tan- 
talizingly invite a 'molecular' interpretation and an analogy with the now well known reso- 
nance V(3872)[0j at the threshold of the meson pair D°D*°. However, unlike the V(3872) 
resonance, which corresponds to an S"-wave meson pair, the suggested resonance would cor- 
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respond to a P-wave charmed meson pair. The appearance of threshold singularities in 
different partial waves can imply a quite nontrivial dynamics of the strong interaction be- 
tween heavy-light mesons. 

In what follows we use the CLEO-c scan data [2] for the e + e~ annihilation cross section 
at six values of the cm. energy (in MeV): 3070, 3090, 4010, 4015, 4030 and 4060. The cross 
section data are available separately for the channels DD, D S D S , DD* and D*D*. Where 
appropriate, these data are for the sum over the pairs of charged and neutral mesons, and 
we assume that the mesons are produced by the isotopically singlet electromagnetic current 
of the charmed quarks, and the isotopic asymmetry, relevant for the channel D*D*, due to 
the isotopic mass difference and due to the Coulomb interaction is taken into account. We 
attempt a fit of the data to a model of a smooth background and a resonance, taking into 
account the interference effects and the effects due to the onset of the D*D* channel at the 
threshold and the channel mixing. In the actual calculation we use, instead of the cross 
sections, the dimensionless production rate coefficients Rk defined as 

a ( e + e - ^ D D) = a (s) 2v 3 D R 1 , a(e + e - -> D S D S ) = a (s) v 3 Ds R 2 , (1) 
a(e + e- -> DD* + D* D) = a (s) 6 R 3 , a(e + e~ — ► D*D*) = a (s) 7 (v'q + v+) P 4 , 

with o" = 7ra 2 /(3s). Here vp, Vd s , v and v + stand for the cm. velocities of each of the 
mesons in respectively the channels DD, D S D S , D*°D*° and D* + D*~ , while for the channel 
DD* + D*D with mesons of unequal mass the velocity factor is replaced by (2p/ y/s) with p 
being the cm. momentum. The extra factors 1, 3 and 7 in respectively the pseudoscalar- 
pseudoscalar, pseudoscalar-vector and the vector-vector channels correspond to the ratio of 
the corresponding production cross section in the simplest model|21 E] of independent quark 
spins, so that the inequality between R 1 , R 3 and _R 4 also illustrates a conspicuous deviation 
from this model. 

The values of the coefficients Ri corresponding to the CLEO-c data at the discussed 
energies are given in the Tabled One can readily see from the table that the rate coefficient 
i?4 is very large near the threshold of the D*D* pairs, a behavior noticed long ago[2], and 
which behavior can make sense only if the production of the vector meson pairs is dominated 
by a resonance, while for the other channels both a resonant and a non-resonant production 
amplitudes are generally present. Accordingly, we parametrize the rate coefficients asjH] 
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Table 1: The dimensionless production rate coefficients Rk denned by Eq.Q calculated from 
the CLEO-c scan data [5]. 



E(MeV) 


3070 


3090 


4010 


4015 


4030 


4060 


Ri 


1.89 ±0.27 


1.18 ±0.23 


0.61 ±0.15 


0.28 ±0.24 


2.92 ±0.27 


3.69 ±0.25 


R2 


34.2 ±8.7 


22.6 ±5.3 


28.9 ±3.2 


23.0 ±5.4 


14.0 ±2.9 


2.7 ± 1.5 


R 3 


49.8 ± 1.4 


45.6 ± 1.2 


45.4 ±0.8 


47.9 ± 1.6 


36.8 ±0.9 


17.5 ±0.6 


R 4 








709 ± 161 


357 ± 14 


81.8 ±3.1 



for k=l, 2, and 3, and 
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where a^, bk, c&, and b are real numbers 1 , and D(E) is the resonant denominator, which 
takes into account the threshold effects due to the strong coupling of the resonance to the 
D*D* channel. The expression for the latter factor at E above the D* + D*~ threshold has 
the form 



D = E — W + - 



_ F , (E-2M ) 3 / 2 (E-2M+) 3 / 2 
1 1 H — r 



(4) 



w 1 ! 2 w 1 / 2 

where Mo (M + ) is the mass of the D*° (D* + ) meson, Wq is the 'nominal' mass of the 
resonance, and T is the width of the resonance decay into all channels except for D*D*. 
The decay rate into the latter final states rapidly changes near the threshold and is accounted 
by the last two terms in the square braces in Eq.Q with w being a parameter with dimension 
of energy, describing the coupling of the discussed resonance to D*D*. At energy E below 
one or both of the D*D* thresholds, the resonance denominator D(E) is found from Eq.p|) 
by analytical continuation. Thus at the energy between the thresholds for the neutral and 
the charged D* mesons, 2 M < E < 2 M + , one finds 



D 1(2M 
4 



- Ef' 2 
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and below both thresholds, i.e. at E < 2 M 



D 



1 

4 



w 



1/2 



Ef'^ + (2_M - Ef' 2 



w 



1/2 



+ E 



w 
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1 The insensitivity of the cross section to the overall phases of the production amplitudes allows one to 
set e.g. the fit parameters and the factor b for the D*D* channel as real. 
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Using the described model of a resonance plus background production of the charmed 
mesons we have made a fit to the data in Table ^ Although the resulting central values 
of the fit parameters (Wo = 4023 MeV, Tq = 67MeV) for the resonance are in a reasonable 
agreement with the data listed in the PDG Tables^ for ^(4040), the statistical significance 
of the fit is quite poor: xV-^dof — 17.6/8. We also have not found any improvement in the 
statistical significance by introducing a linear in the momentum form-factor (1 — p 2 /fJ, 2 ) in 
the amplitudes. Another modification of the described model by introducing a non-resonant 
background for the D* pair production amplitude is not technically feasible right now due 
to lack of data. Indeed, there are in fact two P-wave amplitudes describing the production 
of the D*D* pairs near the threshold: one with the total spin of the meson pair S = 
and the other with S = 2 2 . The resonant production amplitude in the considered model 
corresponds to a fixed composition of these two amplitudes, which composition is generally 
different in the non- resonant part, which thus introduces at least four new parameters. It 
should be however noticed that introducing a non-resonant production of the D* mesons at 
the threshold is neither physically palatable, nor it is likely to improve agreement with the 
data. Indeed, in order to significantly affect the cross section the non-resonant amplitude 
should be comparable with the resonant one, which is very big, as can be seen from the R4 
entries in Table ^ We are not aware of any reasonable mechanism that would explain such 
an enhancement of the non- resonant amplitude. On the other hand the rather large width 
of the -^(4040) resonance implies only a smooth variation of the cross section on the energy 
scale considered here. The data however apparently display a small but sharp 'wiggle' in 
the DD data at 4010 and 4015 MeV and a possible steep rise towards the threshold in the 
D*D* cross section at 4015 MeV. In fact this feature in the data is the main contributor to 
a large y 2 in our fit. Thus introducing a smooth non- resonant term in the D*D* production 
amplitude is very unlikely to improve the situation. 

An obvious effect, which is significantly enhanced in the immediate vicinity of the thresh- 
old for the charged mesons, D* + D*~ , is the Coulomb interaction between the mesons. There- 
fore one can probe this effect as a possible reason for the steeper behavior of the cross sections 
near the threshold indicated by the data. However, we have not found any improvement in 
the statistical significance of the fit by including the Coulomb interaction. In other words 
at the available scan points the Coulomb effect is too weak to be visible in the data. We 

2 It looks reasonable to neglect in the threshold region the third possible amplitude, corresponding to an 
F-wave production with S = 2. 
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nevertheless present here the relevant formulas, which can be helpful if more detailed data 
immediately near the threshold become available. 

Generally, the effect of the Coulomb interaction can be complicated [H] in the resonance 
region by the interference with the strong phase. However due to the relatively large width 
of -0(4040) the variation of the strong phase on the scale of few MeV around the threshold 
(where the Coulomb effect is most essential) can be neglected, and one can neglect such 
interference and estimate the effect of the Coulomb interaction by considering the slow D*^ 
mesons as point charges. In this approximation the Coulomb effect in the D* + D*~ channel 
can be found by adapting the textbook formulas |10j. The corrected expression for R4 takes 
the form 

~i$ + vlF e (E) 



(7) 



\D C {E)\ 2 

and in the expressions for rest of the rate coefficients, as well as in the latter formula, the 
resonant denominator D(E) is replaced by the Coulomb-corrected D C (E), which is described 
as follows. At E > 2 M + , i.e. above the charged meson threshold the term (E-2 M + f/ 2 /w l/2 
in Eq.(jlJ) is replaced as 

(g-2Mjff (£-2M + ) 3 / 2 F ^ 

while below the charged meson threshold the corresponding term in the equations (J5J) and 
(0) is replaced as 

(2M+-E) 3 / 2 {JAU-Ef^- 

where F c (F c ) is the Coulomb factor in the P wave above (below) the threshold. Above 
the threshold this factor is conveniently expressed in terms of the Coulomb parameter A = 



(10) 



M + a/(2k), where k = yM + (E — 2M + ) is the cm. momentum of the D* meson: 

Fc _ 27rA(l + A 2 ) 
1 — exp(— 27r A) 

The analytical continuation of this expression below the threshold can be written in terms 



of the real quantity v = M+a/(2/c) with « = y / M + (2M + - E): 

F c = (1 - z/ 2 ) V (^7T cotTrz/- 2 ln^ . (11) 

The parameter kq in the latter formula can be chosen arbitrarily, and any variation in such 
choice can be absorbed in a shift in Wo and rescaling the resonant amplitude coefficients bk, 
Cfc, and b. 
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The only hypothesis that we have found to be in a statistically significant agreement 
with the data is that the cross section at 4010 and 4015 MeV for the DD channel and at 
4015 MeV for the D*D* channel is contributed by a new resonance. The resonance has to 
be quite narrow, with a width of at most few MeV, in order to not affect the cross section 
at the rest of the scan energies. Within this hypothesis the 4010 and 4015 MeV data for 
these channels are to be removed from the fit, and the remaining data are to be fit within 
the previously described model. At this point we cannot assess from the available data if 
the suggested resonance has any significant coupling to the D*D or D S D S channels. In the 
fit we have made the simplest assumption that such coupling can be neglected and retained 
all the data for the latter channels in the fit. In this way we find that the quality of the fit 
to the described subset of data is quite acceptable: xV-^dof = 3.0/5, with the central values 
of the resonance parameters: Wq = 4019 MeV, T = 65 MeV and w = 9.85 MeV. The plots 
corresponding to this fit are shown in Fig. 1. 




3980 4020 4060 4020 4040 4060 

Figure 1: The plots of the rate coefficients R corresponding to the fit with excluded data 
points at energy 4010 and 4015 MeV for the DD channel and at 4015 MeV for the D*D* 
channel. The excluded points are shown by filled circles. 
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If the deviation of the data points at 4010 and 4015 MeV indeed implies an existence 
of relatively narrow resonance in addition to a broader and more prominent resonance, 
conventionally associated with ^(4040) 3 , an interesting problem arises of explaining the co- 
existence of the two resonances with the same quantum numbers J PC = 1 but apparently 
with strikingly different decay properties. 

One such difference, in the strength of the coupling to e + e~, can in fact be estimated 
from the data. Indeed, if the production of D*D* at E = 4030 MeV is dominated by the 
resonance ?/>(4040), the e + e~ width of the resonance, r ee (4040), can be found from the well 
known formula 

Wga(e + e- -> D*D*) 



r 

1 ee 



3vrr 



D*D* 



(£-iy ) 2 + 7( r o + rW) 2 



where 



1 

D*D* ~ j 



(£-2 M ) 3 / 2 (E-2M+) 3 / 2 



(12) 



(13) 



w l/2 W 1 / 2 

is the (energy-dependent) width of the resonance decay into D*D* pairs. Using the results 
of the fit for the parameters Wo, T , and w, one can readily estimate r Z) *5*(4030) ~ 16 MeV, 
and T ee w 1.7 KeV . 

The coupling to the e + e~ channel of the suggested new resonance X can be very ap- 
proximately estimated from a formula, similar to Eq. ()12j1 applied at the maximum of the 
resonance: 

sj.x) = M ^ (e+ :; - x) - . d4) 

lZ7T 

An estimate of the cross section, associated with the hypothetical new resonance is naturally 

uncertain. The deviation of the data point at 4015 MeV in the D*D* channel is comparable 

to the error, and is about 0.15 of the already small (due to very small phase space) cross 

section cr(e + e~ — > D*D*) pa 0.15 nb. Therefore only about 0.02 of the e + e~ annihilation cross 

section at this energy can be associated with the D*°D*° channel. On the other hand in the 

DD channel the cross section has a minimum at 4015 MeV, apparently due to a destructive 

interference with the background. Thus, if the 'would be cross section' a(e + e~ — >• X — > DD) 

is estimated as the depth of this minimum (i.e. assuming maximal negative interference): 

cr(e + e~ — ■> X — > DD) ~ 0.1 nb, one can arrive at the estimate B ee (X) ~ 10~ 7 . Since, as can 

3 We use the PDG notation ■0(4040) for the dominant resonance, even though its 'nominal' position is 
given at 4019 MeV by the fit. 
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be seen from the plot for R\ in Figure 1, the total width of the X resonance is likely to be 
few MeV, one can estimate the e + e~ width of X in the ballpark of few tenths of eV. 

Clearly, such large ratio of the e + e~ widths (of order 10 4 ) of two closely spaced resonances 
implies that either the mixing between them is finely tuned and results in a great cancellation 
of the coupling of X to e + e~, or the mixing is very small: not larger than 0(0.01) in the 
amplitude. A small mixing can be due to a completely different structure of the two states. 
Indeed the dominant relatively broad resonance can well be a 3 3 Si charmonium state. The 
estimated T ee is close to that predicted for such a state in a charmonium model [7] (T ee = 
1.5 KeV), although the position is somewhat off that in the model (M(3S) = 4.11 GeV). It 
can be also noted that our estimate of T ee as well as the model prediction is somewhat higher 
than the PDG value T ee = (0.75 ± 0.15) KeV for -^(4040). However, given the uncertainties 
and that the data and their analyses are still in flux, we believe that the disagreement 
is not dramatic. In either case it is quite clear that the e + e~ width of ^(4040) is in the 
KeV range, i.e. that typical of a 3 3 Si charmonium state. On the other hand, the new 
hypothetical resonance X can be naturally explained as a D*D* P-wave molecule, similar 
to the S- wave D*°D° molecule X(3872). In particular the small width T ee (X) then agrees 
with the expectedjS] e + e~ width of a P-wave molecule with typical size 0(1 fm). 

The considered here interpretation of the peculiar behavior of the CLEO-c scan data in 
the immediate vicinity of the D*D* threshold in terms of existence of a narrow resonance X 
in addition to ^(4040) can be either confirmed or rejected if more detailed e + e~ data in this 
energy region become available. We would like to point out in this connection that it would 
also be extremely instructive if angular data for the D*D* channel were available. Namely, 
as previously mentioned, the vector meson pairs can be produced in two P-wave states with 
different total spin: S = and S = 2. If the picture of two resonances considered here is 
correct, then it is highly likely that the states X and ^(4040) are coupled to quite different 
combinations of these two spin states. The interplay between the two resonances should 
then result in a rapid variation with energy of a measurable angular correlation, e.g. the 
angular distribution of the pions from the decays D* — > Dn with respect to the direction of 
the initial beams, which distribution is sensitive to both the relative absolute value and the 
phase of the two spin amplitudes [TT]. 

While this work was finalized there appeared the paper presenting the Belle data on 
the cross section of the e + e~ annihilation to the D*D* and D*D final states near and above 
4.0 GeV. This encourages one to expect that more detailed results will become available soon, 
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including possibly the data on the DD annihilation channel. 

This work is supported in part by the DOE grant DE-FG02-94ER40823. 
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